Co-operative density wave and giant spin gap in the quarter-filled zigzag ladder 
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Strong co-operative interactions occur between four different broken symmetries involving charge- 
ordering and bond distortions in the quarter-filled correlated zigzag electron ladder. The ground 
state is singlet, with spin gap several times larger than in the spin-Peierls state of the one-dimensional 
quarter-filled chain with the same parameters. We propose the quarter-filled zigzag electron ladder 
model for several different organic charge-transfer solids with coupled pairs of quasi-one-dimensional 
stacks, the spin-gap transition temperatures in which are unusually high. 
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A characteristic feature of one-dimensional (ID) met- 
als is the Peierls transition, where electron-phonon (c-ph) 
interactions lead to the opening of energy gaps in charge 
and spin degrees of freedom (DOF) and an insulating 
ground state. In the presence of strong electron-electron 
(e-e) repulsive interaction, electronic Hamiltonians de- 
scribing a 1/2-filled band may be reduced to an anti- 
ferromagnetic Heisenberg spin Hamiltonian in which the 
charge DOF are absent. With nonzero spin-phonon cou- 
pling a spin gap (SG) appears again in the so-called spin- 
Peierls (SP) state, in which there occurs lattice dimeriza- 
tion accompanied by the formation of singlet spin-bonded 
pairs of electrons. This mechanism of SG formation is ab- 
sent in the two dimensional (2D) 1/2-filled band, where 
antifcrromagnctism rather than the SP state dominates. 

Away from 1/2-filling, charge DOF are important over 
and above the spin DOF, and the SP transition in ID oc- 
curs only after a metal-insulator (M-I) transition freezes 
the charge DOF within an electronic Hamiltonian. The 
most studied case is that of the 1/4-filled chain, where the 
SP transition is accompanied by charge ordering (CO)i. 
Consider the ID extended Peierls-Hubbard model. 
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tron spin, a and K are the intersite e-ph coupling and 
spring constant, respectively, Ajj+i is the Peierls distor- 
tion of the jth bond, and rij = ^^ c, ^Cj^a is the total 
number of electrons on site j. U and V are on-site and 
nearest neighbor (n.n) Coulomb interactions. Within this 
model, for V < Vc{U) the SP bond modulation is accom- 
panied by 2\iF CO, with site charge densities going as 
0.5 + e , 0.5 + e, 0.5 - e , 0.5 - e . We shall refer to 
this CO pattern as ...1100..., or as a bond-charge density 
wave (BCDW)^ For V > Vc{U), the M-I transition is 
to the CO state ...1010..., and the SP transition involves 
dimerization of the bond distances between the "occu- 
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FIG. 1: (a) The ID ...1100... BCDW SP state. Black (white) 
circles represent large (small) charges. Bond strengths de- 
crease as solid bond > double dotted > single dotted, (b) ID 
1/4-filled SP state with ...1010.. CO. (c) The 1/4-filled rectan- 
gular ladder, (d) BCDW state in the 1/4-filled zigzag ladder. 
Bonds are alternately strong and weak along the stacks, and 
have periodicity 4 along the zigzag interstack direction. 



pied" sites. These two possible CO states are shown 
in Figs. 1(a) and (b). The SP transition temperatures 
in 1/4-filled band quasi-lD organic charge-transfer solids 
(GTS) are typically in the range 10 - 20 K^'S. 

Intermediate between ID and 2D are ladder lattices, 
consisting of two or more coupled ID chains. We will 
distinguish between the two most common ladder geome- 
tries, rectangular (see Fig. 1(c)) and zigzag (see Fig. 1(d)). 
The zigzag ladder is sometimes considered also as a ID 
chain with n.n as well as next-nearest neighbor (n.n.n.) 
couplings. The properties of ladder systems have been 
pursued in great detail in the 1/2-filled electron band or 
spin-model case. SG are found in antiferromagnetic S 
= 1/2 ladders for both even- leg rectangular— and zigzag 
systems^. SGs and superconducting pair-pair correlation 
functions in weakly doped even-leg rectangular ladders 
have been investigated within t-J and Hubbard modelg«. 
Away from 1/2-filling, charge and spin orderings in the 
1/4-filled rectangular electron ladder have been inves- 
tigated within an extended Hubbard Hamiltonian with 
nonzero n.n. Coulomb interaction^. There a zigzag 
CO pattern, in which single electrons occupy opposite 
vertices of n.n. rungs, occurs for sufficiently large n.n 
Coulomb interaction V > V^. SG is possible in this GO 
phase, but with very small magnitude^. 

In the present Letter, we examine the ground state bro- 
ken symmetries in the 1/4-filled band zigzag electron lad- 
der. We show that because of the strong cooperative in- 



teraction between multiple broken symmetries, SGs con- 
siderably larger than those found in the 1/4-filled single 
chain are possible. We then point out that that this 
state is likely realized in certain organic CTS with cou- 
pled chains. 

To understand the BCDW and cooperative effects in 
the 1/4-filled zigzag ladder, we first consider the appro- 
priate one-electron Hamiltonian which includes intersite 
e-ph interactions. 
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In the above, subscripts s and d label respectively in- 
trastack and diagonal zigzag interstack hopping integrals 
and bond distortions (see Fig. 1(d)), and (••) and [••] de- 
note intrastack and interstack n.n. bonds. The broken 
symmetry of Fig. 1 (d) is then simply the 2ki? Peierls dis- 
tortion that is expected when the zigzag ladder is con- 
sidered as a ID chain with n.n. hopping td and n.n.n. 
hopping ts- Unconditional Peierls transition occurs for 
td > 0.5858 tg. We have chosen td/ts = 0.7. Such large 
td have been estimated for the experimental CTS we are 
interested in and are due to the close interchain con- 
tacts involving sulfur atomgSiSiiS. An interesting aspect 
of this Peierls distorted state is that the expected pe- 
riod 4 2]s.p Peierls bond and charge distortions along the 
zigzag interstack direction are accompanied by period 2 
distortions of the same quantities along the stacks. This 
has important ramifications in the presence of e-e inter- 
actions, as we show below. 

We now ask whether this broken symmetry state per- 
sists when Coulomb interactions are added. To the one- 
electron Hamiltonian of Eq. 2 we add an additional term 
containing e-e interactions: 
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U, Vs, Vd are the on-site, intrastack n.n., and interstack 
n.n. Coulomb interactions. We have investigated Hamil- 
tonian H = Hl'^ + _fff| numerically, performing exact 
diagonalization calculations for 16-site periodic zigzag 
clusters and Constrained Path Quantum Monte Carlo 
(CPMC)ii calculations of long (up to 128 sites) open 
as well as periodic zigzag ladders. Investigating the com- 
plete phase space would require varying three Coulomb 
interaction parameters, two hopping integrals as well as 
the different e-ph couplings. We have therefore restricted 
ourselves to the parameter regime appropriate for or- 
ganic CTSi, VIZ., U > 4i„ ts > td, Vd < Vs < U/2. 
In what follows, we have taken ts — 1. The procedure 
for the self-consistent exact diagonalization is the same 
as in Refs. [lfll2l | for ID periodic rings with n.n. hopping. 
For the range of parameters above the ground state we 
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FIG. 2: (a) and (b). Self-consistent order parameters for the 
1/4-filled 16-site periodic zigzag electron ladder (see text). 
Open circles, filled diamonds and open squares correspond to 
f{Bs), f{Bd) and 2An. (c) Normalized An versus U in the 
16-site ID ring with n.n. only hopping (circles), and for the 
zigzag cluster (squares). K = Vd = for [/ = 0, Vs = 1.0 
for all nonzero U in both the ID ring and the zigzag cluster, 
Vd ~ 0.5 in the zigzag cluster, (d) Order parameters f{Bs), 
f{Bd) and 2 An versus Holstein coupling A^. 



find is either completely uniform, or is the broken sym- 
metry state of Fig. 1(d). We evaluate self-consistently 
all (Sf ), {BfA and (rii) and define order parameters (i) 
f{Bs), the absolute value of the difference between con- 
secutive (Bfj) along the stacks, divided by the average 
(Bfj), (ii) f{Bd), the difference between the strongest 
and the weakest (Bf,) along the zigzag direction, divided 
by the average {BfA, and (iii) An, the difference between 
the charge densities of consecutive sites along the stacks. 
Co-operative broken symmetry requires that all three or- 
der parameters are simultaneously nonzero. 

Fig. 2(a) shows the behavior of the three order param- 
eters for fixed U = Q,Vs = l,Vd — 0.5, as a function 
of the dimensionless e-ph coupling As = a^l ^sts, with 
the other e-ph couphng Ad = a'^/Kdts held at 0. All 
three order parameters simultaneously become nonzero 
at As = A^ = 0.045. A^ should be 0+ in the thermody- 
namic limit for the transition to be unconditional, sug- 
gesting that it should decrease with increasing size of 
finite clusters; we verified that A^ is indeed smaller in 16 
site compared to the 8 site zigzag clusters (for 8 sites, 
As = 0.245). Similar behavior are also seen in Fig. 2(b) 
as a function now of Ad, for fixed As — 0.045. Taken 
together. Figs. 2(a) and (b) clearly indicate the co- 
operative nature of the transition to the BCDW state 
in the 1/4-filled zigzag ladder. 

The broken symmetries within the correlated zigzag 
ladder are far stronger than in ID. To demonstrate this 
we compare An in the two cases as a function of U with 



e-ph coupling A = a^/Kt = 1.28 for ID and A^ = 0.125, 
Xci — Q for the zigzag ladder. The distortion amplitudes 
decrease rapidly with U in ID (see below), which neces- 
sitates the larger A in this case. In Fig. 2(c) we have 
plotted the An normalized by its value for t/ = for 
both cases. For all nonzero U we have chosen V = 1 for 
the ID ring and Vg ~ I, Vd — 0.5 for the zigzag clus- 
ter. The distortion in the zigzag ladder is affected very 
weakly by e-e interactions because of the co-operative 
nature of the transition here. As pointed out above, the 
charge and bond periodicities along the stack direction 
in Fig. 1(d) are Akp, even as these have 2ki;' periodicities 
along the zigzag interstack direction. In purely ID the 
Akp distortions get stronger with e-e interactions, while 
the 2kF distortions get weaker. There is thus a tendency 
to cancellation of these effects in the zigzag cluster, with 
the distortion amplitude remaining the same. 

As one goal of our work is to apply our theory to or- 
ganic CTS with molecular sites we have also investigated 
the effect of Holstein electron-molecular vibration cou- 
pling: 
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Here /3 is the intrasite e-ph coupling with corresponding 
spring constant Kf^, and Vi is the amplitude of the in- 
ternal molecular vibration. In Fig. 2(d) we have plotted 
the order parameters for H = Hll + i/|| -f Hp against 
A^ = 13^/Kfits , with fixed \s ^ Xd = 0.045. AU three 
order parameters again increase with Xp. Our conclu- 
sions regarding cooperative broken symmetries remain 
the same whether or not the Holstein interaction is in- 
cluded. 

We have verified the coexisting broken symmetries for 
larger lattices by performing CPMC calculations for long 
zigzag clusters with open boundary condition (OBC). In 
ID, the central region of long open chains exhibit spon- 
taneous charge or bond-order distortions, even with uni- 
form hopping integrals and zero e-ph coupling e^'^ . We 
have performed CPMC calculations for a 64-site zigzag 
open ladder with Xs = Xd = Xp = Q, and for the same 
[/, V;, t, and td as in Fig. 2(a). With OBC, f{B,)J{Bd) 
and An depend on the locations of the sites being con- 
sidered, and we have therefore plotted the charge densi- 
ties (rii) and the bond orders {Bfj) at chain centers in 
Fig. 3. Simultaneous A]<ip dimerizations of charge and 
bond along the stacks (Figs. 3(a) and (b)) and simul- 
taneous 2kj7 modulations of these quantities along the 
diagonal zigzag direction (Figs. 3(c) and (d)) are seen, in 
complete agreement with the results of Fig. 2. 

The SG in the distorted zigzag 1/4-filled ladder is due 
to the formation of a strong 1-1 local singlet bond in the 
...1100... zigzag BCDW. We have calculated the SG in 
the BCDW state of the zigzag ladder using CPMC. Our 
calculations are for long periodic zigzag clusters, for the 
same parameters as in Fig. 3. For comparison, we also 
evaluate the SGs in the ...1100... BCDW state of a ID 
periodic ring (Eq. (1), Fig. 1(a)), with U — 6, V — 1, t 
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FIG. 3: Charge densities and bond orders in an open 64-site 
1/4-filled zigzag ladder, as computed by CPMC. Site indices 
i are arbitrary and correspond to the central region of the 
open ladder, along the stacks in (a) and (b), and along the 
interstack zigzag bonds in (c) and (d). 



= 1, using the Stochastic Series Expansion Monte Carlo 
method^**. The reason for choosing the ...1100... state 
over the ...1010... CDW-SP state is that the SG in the 
former, with shorter 1-1 bond, is larger. Finite periodic 
clusters do not have spontaneous distortions while self- 
consistent CPMC calculations are not possible for the 
large clusters we need to investigate to obtain the SGs in 
the thermodynamic limit. Thus the BCDWs can be gen- 
erated only by externally imposing bond or charge dis- 
tortions. The co-operative nature of the zigzag BCDW 
demonstrated in Figs. 2(a) imply that either of the two 
bond distortions or the CO is suitable for our purpose, 
as the other two broken symmetries will be generated 
spontaneously. Since the modulations of the hopping in- 
tegrals are different in the ID chain and the zigzag ladder, 
we choose to externally impose COs that are congruent 
with Figs. 1(a) and (d). This is achieved by adding a site 
energy component J2i ^i'^i to the As = Ad = limit of 
the Hamiltonian in Eq J^Jl. The site energies added are 
+ |e|, -|-|e|, — |e|, — |e| along the zigzag bonds as in Fig. 1(d) 
and along the linear chain as in Fig. 1(a). We evaluate the 
singlet-triplet gaps A^ using CPMC for periodic zigzag 
ladders with N = 32, 64 and 128 sites, and for periodic 
ID rings with N = 32, 64 and 96 sites, for several different 
e. For each e, the Ao- and An at A'^ — > oo are found from 
extrapolations against 1/A^. In Fig. 4 we have plotted 
the extrapolated A^. against An for both the ID system 
and the zigzag electron ladder. The finite size scaling of 
Acr for the zigzag ladder for one value of e is shown in 
the inset. For the same An, Aa- is several times larger in 
the zigzag electron ladder than in ID (Aa- can be plot- 
ted also against |e|; the same large difference between the 
zigzag ladder and ID is obtained.) The large SG in the 
zigzag ladder is a direct consequence of the strong local 
interstack 1-1 singlet bond, which in turn results from 
stronger distortions than in ID. 

Several recently discovered 1/4-filled band CTS are 
very likely structurally zigzag ladders. These systems 
consist of pairs of ID stacks of organic donor moleculesi^. 
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FIG. 4: Finite-size scaled spin gap Ao- versus An in the 1/4- 
filled band zigzag ladder (filled symbols) and in the ID chain 
(open symbols). Lines are guides to the eye. Inset shows the 
finite size scaling of Ao- in the zigzag ladder for |e| = 0.3. 



with strong intrapair interchain couplings, and weak in- 
terpair couplings ^^'^^. They undergo M-I transition at 
Tm-i > 200 K, followed by an insulator-insulator (I- 
I) transition which is accompanied by the opening of 
a SG at lower temperature Tso^'^'^^^- The peculiar- 
ity of these paired-stack systemsSi^iiSiiS is that while 
Tgp - 10 - 20 K in the ID CTSSi^, Tsg -^ 70 K in 
(DT-TTF)2Au(mnt)2^ and is even higher at 170 K in 
(BDTFP)2PF6(PhCl)o.3ia. The very large Tsg, and 
therefore the large SG in these systems are highly un- 
usual. 

To explain the large SG in these systems, an ef- 
fective two-leg rectangular spin-ladder model has been 
proposedSiSiiSiiS. It is known experimentally that each 
individual 1/4- filled stack dimerizes below Tmi^- Within 



the dimerized rectangular ladder lattice motif, if each 
dimer unit along the legs containing one localized spin 
is assumed to be equivalent to an effective single site 
(see Fig. 6 in Ref. |3), then an 1/2- filled rectangular 
spin ladder is obtained. For strong enough interstack 
spin coupling, large SG is now obtained. The mapping 
from the 1/4- filled dimerized rectangular electron ladder 
to the 1/2-filled spin ladder, however, has not been for- 
mally proved. 

We believe, however, that the zigzag electron ladder is 
the more appropriate model for these systems for a num- 
ber of reasons. The crystal structures of the materials 
(see Fig. 2(a) in Ref. Hand Fig. 2 in Ref. [13) indicate that 
each molecular site is coupled to two molecules on the 
partner stack, as would occur in the zigzag ladder. Quan- 
tum chemical calculations of hopping integralsi^ support 
this viewpoint. Most importantly, it has been suggested 
from EPR linewidth studies that transition to a CO state 
might be occurring in (BDTFP)2PF6(PhCl)o.5^°. This 
last result, if correct, will be in agreement with the zigzag 
electron ladder model. The key difference between spin- 
ladders and 1/4-filled electron-ladders is that CO is ab- 
sent in the former, while it is a prerequisite to large SG 
in the latter. Experiments have previously demonstrated 
CO in inia as well as 20^° GTS, while very recently, the 
occurrence of the ...1100... BCDW has been experimen- 
tally confirmed in a CTS^. We predict that experimental 
investigation of the coupled-stack GTS will find evidence 
for CO. 
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